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Abstract 

The second-order hydrodynamical description of a homogeneous conformal plasma that undergoes a boost- 
invariant expansion is given by a single nonlinear ordinary differential equation, whose resurgent asymptotic 
properties we study, developing further the recent work of Heller and Spalinski [Phys. Rev. Lett. 115 , 072501 
(2015)]. Resurgence clearly identifies the non-hydrodynamic modes that are exponentially suppressed at late 
times, analogous to the quasi-normal-modes in gravitational language, organizing these modes in terms of 
a trans-series expansion. These modes are analogs of instantons in semi-classical expansions, where the 
damping rate plays the role of the instanton action. We show that this system displays the generic features 
of resurgence, with explicit quantitative relations between the fluctuations about different orders of these 
non-hydrodynamic modes. The imaginary part of the trans-series parameter is identified with the Stokes 
constant, and the real part with the freedom associated with initial conditions. 
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I. INTRODUCTION 

Resurgent asymptotics using trans-series is a powerful way to extract physical information 
from asymptotic expansions, systematically including exponentially small contributions that are 
typically neglected in traditional asymptotic analysis a-la Poincare [1]. This formalism is well 
developed for ordinary differential equations, both linear and nonlinear, and to a lesser extent 
some results are known for partial differential equations. Physical applications have included 
problems in fluid dynamics [2], exact WKB [3-5], matrix models and strings [6, 7], and quantum 
field theory [8-10]. There is also a well-developed literature concerning the Painleve transcendents, 
which themselves have many physical applications, being the nonlinear analogues of the familiar 
special functions of linear physics [11]. A common thread is the goal to incorporate, in a controlled 
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numerical and analytic manner, exponentially small corrections into asymptotic expansions, in 
such a way that the trans-series encodes its proper analytic continuation properties [12-14]. 

In this paper we discuss a new example in physics, in the context of hydrodynamics. This is 
motivated by an interesting recent paper by Heller and Spalinski concerning resummation of the 
gradient expansion in conformal hydrodynamics [15]. This example is relevant for hydrodynamic 
studies of heavy ion collisions [16, 17], and also addresses fundamental issues of the nature of 
the hydrodynamics expansion. High orders of the gradient expansion of the linearized hydrody¬ 
namic equations have been studied in [18, 19]. Here we investigate the gradient expansion in the 
full nonlinear system. Our choice of conformal hydrodynamics, and the boost invariant regime 
thereof, allows us to reduce the hydrodynamic equations from a set of coupled nonlinear partial 
differential equations to a single nonlinear ordinary differential equation. We study in pedagogical 
detail the asymptotic properties of this nonlinear equation, using resurgent trans-series. We show 
that this equation displays resurgent trans-series relations between different non-perturbative sec¬ 
tors, and that the asymptotic hydrodynamic expansion encodes detailed quantitative information 
about non-hydrodynamic modes. These results confirm the relevance of the resurgence approach. 
The ultimate goal is to extend these lessons from the reduced one-dimensional system to the full 
hydrodynamical system away from the boost-invariant and/or conformal limits, which are presum¬ 
ably non-integrable, and also to the corresponding nonlinear partial differential equations arising 
through the AdS/CFT correspondence [20-26]. 

II. NONLINEAR EQUATION FOR BJORKEN FLOW IN CONFORMAL HYDRODY¬ 
NAMICS 

In this section we briefly recall the relevant notation for describing second order hydrodynamics 
and the steps to reduce it to a single nonlinear equation in the restricted case of Bjorken flow in 
conformal hydrodynamics. The general second order dissipative terms for relativistic, conformal 
hydrodynamics have been derived in [27]. The equations of hydrodynamics stem from the conser¬ 
vation of the energy-momentum tensor^, = 0. In order to describe the energy-momentum 

flow we use as our hydrodynamic fields the energy density in the local rest frame of the fluid, £, and 
the fluid four-velocity, with u^Uf^ = — 1. The covariant expression for the energy-momentum 

^ We assume the absence of other conserved charges. 
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tensor in terms of these hydrodynamic fields is 


=p{g)g^^’^+ {p{g)+S)u^^u’^+ ( 1 ) 

where p(g) is the pressure (related to the energy density by the equation of state), and is the 
dissipative part that includes the viscous corrections. Conformal symmetry implies = 0, which 
determines the equation of the state to be p{g) = g/{d — 1), for d space-time dimensions. The 
dissipative part, II^^, is symmetric, transverse (i.e. = 0), and for conformal fluids, traceless. 

In the hydrodynamic limit where one considers only the long wavelength, small momentum modes, 
the dissipative corrections that constitute 11^^ are given as an expansion in space-time gradients 
of the hydrodynamic fields g and [27-29]. This derivative expansion contains all terms, at a 
given order, which are allowed by the nnderlying symmetries of the fluid. At first order in the 
derivative expansion, there is only one term allowed by Lorentz and conformal symmetries, the 
shear viscosity rj{g): 

C* order : = -t] = -2r] < > (2) 


Here we use the standard notation [15, 27, 29] to denote the symmetric, transverse projection of a 
tensor 




( 3 ) 


where -I- is the transverse projection operator. 

At second order, there are five new transport coefficients [27]. One characterizes the relaxation 
of the energy-momentum tensor, one characterizes the conpling to space-time curvatnre, and the 
other three describe non-linear couplings to the fluid velocity. Among these latter three, two are 
related to couplings to vorticity. In this paper we study homogeneous and boost invariant flow, in 
which case we can neglect the terms coupling to curvature and vorticity, so only two of the five 
possible terms contribute: 


2”'^ order 


+ 77Tn ( ^ 


d- 1 


+ Aicr|' 




+ ... 


( 4 ) 


The relevant second order transport coefficients are rn and Ai. Using conformal invariance we 
can parametrize the dimensionful quantites in units of local temperatnre T{x). For example, the 
energy density scales as T oc T^. For the transport coefficients, we adopt the parameterization 
used in [15] 

rj = CnS , rn = ^ , Ai = Cx^ (5) 


where s oc ^ is the local entropy density, and Cr, C\ and Cr^ are dimensionless numbers. 
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A. Bjorken flow 


We consider boost invariant flow, also known as Bjorken flow [30]. It describes a homogeneous 
fluid expanding longitudinally along a fixed direction, say z. The original motivation of Bjorken 
flow was to capture the essential physics of the space-time evolution of the quark-gluon plasma 
produced in heavy ion collisions, which is modeled by a fluid expanding between the highly Lorentz 
contracted “sheets” of nuclei, moving away from each other almost at the speed of light, shortly 
after the collision. With boost invariant initial conditions the expansion is boost invariant, in 
other words the system looks the same in all inertial frames. It is convenient to work with proper 
time, r, and rapidity, C as the space-time coordinates. They are related with the usual Minkowski 
coordinates z and t as 


t = T cosh ( 


z = T sinh C 4^ 


T = \/t^ — z'^ 


C = tanh ^{t/z). 


( 6 ) 


The metric is ds"^ = —+ dx\ and the non-vanishing Christoffel symbols are T^^ = r“^, 
and = r. In these coordinates the fluid velocity is constant, = 1, = 0, and 

the conservation equation, = 0, along with Eq. (1) and the conformal equation of state, 

p = £/{d — 1), reduce to 


d£ d 


T-$ = 0 


( 7 ) 


where [27]. At leading order where the viscous terms $ are neglected, the energy density 

has the asymptotic large r behavior 


T ~ r '*-1 -|- ... , ~ T . , d = 4^ 


( 8 ) 


Here the ellipses denote the viscous corrections. Equation (8) is the leading order result that 
describes how the system cools as it expands. The viscous corrections can be calculated by inserting 
(4) into (7) and solving as an expansion at large proper-time r. We instead follow a Miiller-Israel- 
Stewart-like approach and promote the dissipative part of the energy-momentum tensor H^'^ to 
an independent hydrodynamic field.^ Using the first order relation (2) we can rewrite (4), also to 
second order, as 




rjTu 


(9) 


^ The Muller-Israel-Stewart approach [31] was designed to obtain a set of hyperbolic differential equations which 
have causal solutions as opposed to the hydrodynamical equations which are not hyperbolic and have acausal 
propagation of certain modes. These modes are beyond hydrodynamics, but eliminating them has advantages in 
numerical simulations. 
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and obtain a relaxation equation for with relaxation time rn- This procedure generates an 
all-orders derivative expansion which agrees with the ordinary hydrodynamic expansion up to third 
order, beyond which more transport coefficients should be taken into account. 

In the homogeneous, boost invariant limit the relaxation equation (9) further reduces to a 
non-linear equation 


^ , f-, , dm 


, , d — 3\ Ai,0 

-2 




d — 1 ) T 


( 10 ) 


where the transport coefficients rj, m and Ai are functions of the energy density £. The technical 
problem now is to solve the two coupled nonlinear equations (7) and (10). 

To illustrate the ideas of resurgence in hydrodynamics, we consider as our central object the 
all-orders expansion generated by this Miiller-Israel-Stewart-like treatment of the second order, 
conformal, boost invariant hydrodynamics. In d = 4 the resummation of this expansion is studied 
by Heller and Spalinski in [15]. We show that the usual derivative expansion encodes much more 
information, such as the non-hydrodynamical modes, than one might naively anticipate, a result 
which we argue has implications for more complicated and phenomenologically relevant models as 
well, such as those relating hydrodynamics to gravitational systems [21, 22, 24, 26]. 


B. Scale invariance and the nonlinear hydrodynamic equation 


Using the parametrization of the transport coefficients (5) dictated by conformal invariance, 
equations (7) and (10) combine into a single highly nonlinear equation for the temperature T(r): 






T 

1 


d- 1 


{d-2)Cn T ' r2 

, (d-3)CAA dCr-2(d-2)C„ ^ 

^ (d-2)cj '' ^ 


d- 1 


( 2id-3)C,\_,. 

( 11 ) 


(d-l)2 

where T = dT/dr. As a consequence of the underlying scale invariance of the system, Eq. (11) 
is invariant under rescaling r —)• ar, T —)• a~^T. Therefore we can integrate Eq. (11) to obtain 
a first order equation. Scale invariance also suggests introducing a local proper-time variable, w, 
measured in units of local temperature (in mathematical terms, w is known as the “Ecalle time”): 


w = tT(t) 


( 12 ) 


6 













and an associated dimensionless characterization of the temperature, through a function f{w) 
defined as^ 




dlogT 


(13) 


dlogr 

With these definitions, Eq. (11) reduces to a single first order nonlinear equation for the dimension¬ 
less variable f{w) as a function of w, the proper time measured in the units of local temperature. 

CrWf{w)f{w) + CrWf'{w) + ^/(^) + dCrf{w) + ^ 


, I , {d-3)Cx\ , 2dCr,, , , dCr 

I'+ (d3^J +— 


- 2{d-2)Cn 


= 0 


(14) 


{d-iy 

(Here we study this equation in d = 4, but we note that d = 3 also has novel features.) For d = 4, 
we reach the central equation to be studied in this paper: 


CrWf{w)f'{w) + Crwf{w) + (^^ + wf{w) + iCrfiw) + ^^^/^(^^) 

+ 3 ^ ° 

Up to the simple shift, /ours = (/theirs — 1) noted above, (15) is the equation analyzed in [15]. 

It is clear that after the rescaling w —>• CrW, Eq. (15) only depends on the ratios CrjfCx, and 
Cr^jCr- Physically, this is simply because Cr characterizes the relaxation time, using w now as the 
time variable, and so is naturally absorbed into w. However, to facilitate direct comparison with 
previous results [15], we will not adopt this rescaling. And in what follows, for concrete numerical 
illustrations, we will also use the parameters that are associated with AA = 4 super Yang-Mills 
theory [15, 27, 29]: 




& = 


2 - log 2 




( 16 ) 


III. FORMAL LATE TIME EXPANSIONS 


In this section we analyze the hydrodynamic expansion generated by the nonlinear equation 
Eq. (15). The hydrodynamic regime is identified with large proper time, or equivalently large w. 
Mathematically speaking the large w expansion is a formal series which is divergent, asymptotic 
and non-Borel summable. Such asymptotic behavior is characteristic of gradient expansions of 
effective actions [32]. This is by no means a bad thing: the asymptotic nature of the expansion 
actually encodes important physical information. 

® Note that our / differs from the / in [15] by: /ours = /theirs — 1. 
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A. Late time hydrodynamic expansions 


The late proper-time hydrodynamic expansion"^ of f{w) is an expansion in inverse powers of w. 
Since w ~ in d = 4, this expansion translates into a late proper-time expansion for T(r) using 
(12, 13). We start with the formal series ansatz 




( 0 ) 


w 


(17) 


k=0 


The meaning of the superscript (0) will become clear in Section IV, when we include the terms 
beyond hydrodynamics. With this ansatz (17), the nonlinear equation (15) generates a recursion 
relation for the the coefficients 

k / ^ \ k—1 

do) 

'k' 


Cr t. (*=' -1) VT /1°1. -Cr(k-i) + (i + ((^) /f+4a V T-V-i/i 

fc'=0 ^ k'=0 

E T-7/f + ^ (‘ + ^) = “ 08) 


These can be solved iteratively. The equation that determines is quadratic. Therefore there 
are actually two solutions to the recursion relations, leading to two different large w expansions: 


f-\w) 


1 4C„ 1 

-h ^ 

3 


8Crj{Cx-Cr) 1 , 16C'^(2(CA-a)"-3C^a) 1 


1 


-X 1 + 


9 w 

2 a 


27 




+ 


81 


W'^ 


Cx 


4a 


1 - 


4a a \ 1 


''rj 


-h • • • 

W 


+ ... (19) 

( 20 ) 


In terms of physical quantities, such as the proper time r and the local temperature T(r), it is the 
former expansion, f^\ which leads to the familiar hydrodynamic expansion T(r) ~ -|- 

.... On the other hand, the latter expansion /_ leads to T{t) ~ r ® V , which does 

not connect to the ideal hydrodynamic result at late times. Furthermore, the expansion f^\'w) 
in (20) is unstable, in a sense explained below, so we concentrate on the expansion from now 
on. With this understood, for notational simplicity we drop the subscript “ , and simply write 

j(o)(-u;) for the function with late time expansion in (19). 


As terminology, we refer to the all-orders expansion (17) generated by the Miiller-Israel-Stewart-type analysis 
described in Section II as the “hydrodynamic expansion”. We will see that a consistent asymptotic analysis of the 
basic nonlinear equation (15) requires the addition of further terms, which go beyond hydrodynamics. 












B. Divergence, Borel ambiguities, and exponentially suppressed terms 


It is straightforward to show that the expansion (19) is divergent. The late terms in the 
expansion grow factorially fast, and their leading large-order growth can be characterized as 




( 0 ) 


r{k + i3) 


k 


oo 


( 21 ) 


'k gk+P 

for some real numbers S and /3, whose physical meaning will become clear shortly. Furthermore, 

with the J\f = 4: parameters in (16), the parameter S turns out to be positive. There are three 

important and interrelated points related to this divergent behavior [1, 33]. First, the hydrodynamic 

expansion (17) is a typical asymptotic expansion: for a hxed large value of w, the terms at 

first decrease for increasing low values of fc, but eventually start growing beyond some value of k, 

say k*. This transition happens roughly when d{fj^^w~^)/dk = 0, or k* ~ Sw. One can obtain 

k* 

exponential precision by the procedure of “least-term-truncation”, truncating the series at fc ~ A:*. 
The associated inherent error is exponentially small 




^ e 




Second, one could try to apply Borel summation: 




“ F(n + /3) _ /■“ du 


71=0 


^n+/3 


W 


E 

n=0 


W 


^-Swu ^n+/3 


f 


du 


~ ; — e 

u 


—Swu 


1 — U 


( 22 ) 


(23) 


However, when S' > 0, as is the case here, the Borel integral has a pole at u = 1, leading to an 
ambiguous imaginary part 




(24) 


Notice that this is the same order of magnitude as the least-term-truncation error in (22). The 
ambiguity (24) is due to the ambiguity of deforming the integration contour around the pole. The 
existence of such an ambiguous term is problematic for several reasons: it has an undetermined 
sign, and furthermore it is pure imaginary. However the function f{w) should be real, since it is 
related to the physical temperature. It is clear that in order to fix this ambiguity we need to add 
to (w) an exponentially suppressed term of the same order. 

Third, one can consider linearized perturbations around the solution This is achieved 

by an ansatz —>• + df^^\w). Inserting this ansatz into Eq. (15), keeping the first 

two terms in the w~^ expansion of f^^\w) in (19) and the linear terms in Sf^^\w) leads to the 
equation 


(25) 
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which has the solution^ 


/S = ( 26 ) 

Notice that the linearized perturbation in (26) has the same functional form as the least-term- 
truncation error (22), and the ambiguity (24) of naive Borel summation. This is not a coincidence, 
as all three phenomena are manifestations of the asymptotic character of the hydrodynamic ex¬ 
pansion (17). In fact, by general arguments [1, 33] the constants S and f3 appearing in the leading 
large order growth (21), and hence also in the error (22) and the ambiguity (24), are exactly the 
same constants that appear in the linearized perturbation (26). Thus we deduce that 


5/(0) (, 


w) 


exp 


2a 


■ w 


= J_ . _ - 2Ca 

2Cr ’ ^ Cr 


(27) 


which becomes now a numerical prediction for the leading large-order growth (21) of the coefficients 
of the hydrodynamic expansion, which are generated from the recursion relations (18). Numerically, 
for the A7 = 4 SYM plasma parameters in (16): 


S 


3tt 

2 - log 2 


7.21181 


/3 = 


3 

4 - 2 log 2 


-1.1478 


(28) 


In Section V we confirm the consistency of these arguments with great precision. 

The conclusion is that by itself the hydrodynamic series expansion in (17) is merely a 

formal expression. In order to promote it to a well-defined function it is necessary to enhance it 
by adding some exponential terms that go beyond the hydrodynamic expansion. These terms are 
not arbitrary: they are highly constrained. The linearized perturbation is only the first element 
of an infinite set of exponential corrections. These exponential terms are required to render an 
unambiguous answer for the actual solution, f{w), of the differential equation. The full expansion 
that contains all these exponential terms in addition to the formal series expansion (17) is known 
as a “trans-series” [1]. In the next section we demonstrate how such a trans-series is constructed. 


IV. TRANS-SERIES EXPANSIONS 

The leading exponential correction to the perturbative series, is only the tip of the 

iceberg. Just like the hydrodynamic series itself, the first exponential correction is also 

® The same argument for the other expansion f^\ in Eq. (20), leads to 5f^\w) ~ exp [^~ 2 c (C^-Cx) 

2([/3_4C^C' 

B = — Y T ^ ^ .I’ ^ . Notice that for the Af — 4 parameters (16), the exponent is positive, which is related to 
the unstable nature of the expansion in (20). 
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an asymptotic series: the factor e~^'^ is multiplied by another formal series in w~^, which 
we denote as and which is also divergent, and with the same leading rate of growth 

(21) as in the hydrodynamic series. By the same arguments, we conclude that there must be a 
correction of the form Furthermore, this second exponential term is also multiplied 

by another asymptotic series, denoted by which requires the existence of corrections of 

the from This pattern continues ad infinitum. (This is also clear generically from the 

nonlinearity of the underlying differential equation.) In order to obtain an unambiguious answer, 
all these exponential terms, and their fluctnations, must be included in the answer. The resnlting 
combined object is known as a “trans-series”More importantly, the theory of resurgence predicts 
that these various asymptotic series, which make up the full trans-series, are related to one another 
in extremely intricate ways. These inter-relations will be demonstrated explicitly in Section V. 

We are led to the following trans-series ansatz generalization of the formal hydrodynamic ex¬ 
pansion (17); 

f{w) ~ + + + (29) 

OO 

~ (30) 

n=Q 

This is a sum over powers of an exponential factor (the analog of the “instanton fugacity” in semi- 
classical expansions in quantum mechanics or quantum field theory, with S being the instanton 
action) 

Ciw) = wf^e-^^ (31) 

each multiplied by a formal large w series (the analog of the pertrnrbative fluctuations about the 
instanton sector): 



In the trans-series expansion (29, 30), a is the trans-series expansion parameter, which is in general 
a complex number. As discussed below, its imaginary part will be fixed by resurgent cancellations 
associated with the reality of the trans-series, while its real part is a free parameter related to the 
initial conditions of the ODE. At this stage, it simply connts the “non-perturbative” order of the 
trans-series expansion. 

® More general trans-series also include powers of logarithms, and possibly iterations of powers, exponentials and 
logarithms [1, 13, 14]. 
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Inserting the trans-series ansatz (29) into the nonlinear ODE (15), we obtain a set of recursion 
relations by equating to zero the coefficient of each ^ > q n > 0. Setting 

n = 0 leads to the recursion relations (18) that generate the hydrodynamic series (19) studied in 
the previous section. Setting n = 1, we obtain the recursion relations that generate the asymptotic 
series {w) that multiplies the first exponentially suppressed term, or the first non-hydrodynamic 
series. We refer to the asymptotic series as the non-hydrodynamic series”. 

-Cr P [/Ih, (t' - 1) + fZ. (S/<;> + {k'-fi-l) /“ .)] + (l + )|) fi" 


-Cr{Sf<'' + (k-iS-l)f<Z) +4C. 


fc-1 


r f(i) f(o) , ^(0) .(1) 

[jk-k'-l-Zk' Jk-k'-lJk' 




3Ca 

2C, 


K 

E[/. 


( 1 ) f(o) , AO) Ai) 
k-k'Jk' '^Jk-k'Jk' 


V i.,_ 




fc'=0 

SCrAl) 


A > - 0 
Jk-i — ^ ■ 


(33) 


k'=0 


The first two terms, k = 0 and k = 1, in these recursion relation determine the constants S and 
/3 to be as in (26). The leading coefficient is a free parameter, and the remaining coefficients, 
[k > 1), are determined by Eq. (33), and are all proportional to Therefore, without loss 
of generality we can normalize = 1, and characterize the freedom in terms of the constant a. 
This normalization fixes all the other coefficients uniquely. The leading non-hydrodynamic 


series is 


. 2{C^^-Cr,iC^-6Cr) + 2CxiCr-C^)) 

p >{w) ~ 1 -—- 

3Ct w 


+ ... 


(34) 


The higher hydrodynamic series are determined by similar recursion relations. For completeness 
we give the full recursion relation for arbitrary n and k: 


-c. EE// 


(n—n') 


k-k' 


n'=0 fc'=0 


n' .S + {k' - n'/3 - 1) fpPPj - a (n.S /f ^ + (k - n/3 - 1) f^J^j 




Cr^ 


n fc—1 


(n—n') f{n') 3 C\ 


+ 


n'=0 k'=0 


2 a 


EE/ 


(n—n') p(n') 


k—k' 


/' ^ + 


8 a 


An) 

Jk-1 


^ n'=0 k'=0 


+ 2 + 2 ^ ^ ^k,0^n,0 + g {Ct — Cri) ^k-l,0^n,0 — 0 (35) 

Once the free parameter /q^\ or a, is fixed there is no more freedom, and all the coefficients 

are determined completely by the recursion relations (35). This is perhaps surprising, but it is 

generic [1]. The dependence of the trans-series on a is exactly as written in Eq. (30), i.e. cr enters 

the trans-series as For example, the second non-hydrodynamic expansion is 

^.2), ^ 3(Ca - Cy) -2C3 + C^{ACx - llCr) + 2Cy (C| + ACxCr + 2CPj + 2C|(a - 2 Ca) 
/ W) ~ - AA -^-- 


2G 


CyCrW 


+ ... 


(36) 
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and the third non-hydrodynamic series is 

/»>(«,) ~ ~ ~ + 4C^ 

4C,C?{4Cj - 5C,) + 3(C, - C;i)(a, + Ca)(C, - 2Ca)( 3C, - 2Cj)] + .., (37) 


The remaining coefficients can be generated in a straightforward fashion, but are rather cumber¬ 
some to write. 

More physically, these non-hydrodynamical series correspond to exponentially damped modes 
multiplied by a gradient expansion due to viscous terms. Recalling w ~ these modes con¬ 
tribute to expansions of physical quantities such as the local temperature, or energy density as 


^(w) 


2n/3 _„_3_^2/3 ^ — 

q- 3 Q 2 Ct ~ 7" 3 g 2 To 


(38) 


Note that for the AA = 4 parameters (16) these modes introduce a transcendental power (i.e. /?) of 
T in the gradient expansion. It is illustrative to compare the non-hydrodynamic modes to instanton 
contributions in quantum mechanics and quantum field theory. The coefficient S that controls the 
damping is analogous to the two-instanton action in the dilute instanton gas picture, where the 
ambiguities that arise from the divergence of perturbation theory are cured by non-perturbative 
corrections from higher (multi-)instanton sectors [4, 5, 9]. Here the non-hydrodynamic modes play 
the role of instantons. 


V. LARGE ORDER BEHAVIOR, STOKES CONSTANTS AND BOREL ANALYSIS 

The reality condition on the trans-series f{w) in (29, 30) means that cancellations must occur of 
imaginary terms generated by Borel summation of the different (w) expansions. As described 
in the comprehensive analysis of real trans-series by Aniceto and Schiappa [13], this leads to an 
infinite hierarchy of relations between the expansion coefficients of different n sectors. In 
this Section we analyze these relations numerically, and confirm that these resurgence relations 
do indeed hold. This also allows us to deduce the Borel plane structure of the associated Borel 
transforms, which reveals some interesting branch-cut structure. 


A. Hydrodynamic derivative expansion: 


Resurgence predicts [13] that the leading large-order growth of the expansion coefficients 
in (21) receives subleading corrections of the form: 

5 


( 0 ) 


AO) 

Jk 


r{k + f3) / ( 1 ) 

‘^'27rz5^+/3 + 


k + P-1 


+ 


{k + /3-l)ik + ^-2) 




+ ... . (39) 
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This is an example of the resurgent relations between two sectors of the trans-series, namely 
between the hydrodynamic and the first non hydrodynamic series Recall that the 

constants S and /3 appearing in (39) are known, being given by (27). The subleading terms are also 
known, as they are completely determined by the low order coefficients, /q^\ f 2 ^\ ..., of the 

first non-hydrodynamic series in (34). It is therefore possible to determine numerically the 

overall constant Si, a “Stokes constant”^, by generating many coefficients from the recursion 
relations (18). It is simple to generate many thousands of these coefficients, which provides great 
precision. 

In Figure 1 we plot the ratio of the large order expression (39) to the actual coefficients. This 
clearly demonstrates the remarkable precision of the resurgence relation (39). Indeed, as shown in 
the second plot, when all three terms on the right-hand-side of (39) are included, the agreement 
with the prediction of the large-order growth is at the one percent level already by the 10**^ term in 
the expansion. From direct comparison of the large order expression (39) with the exact coefficients 


(large order) ♦ 1-term 



(large order) « 3_terms 

(exact) 
l.lOf 

1.05 

1.00 . . 

0.95 

0.90 

0.85 

2 4 6 8 10 12 14 


FIG. 1: The first plot shows plot the ratio of the large order expression (39) to the exact coefficients 
as a function of k. The blue, green and red points refer to the inclusion of the zeroth, first, and second 
subleading terms in (39). The second plot shows a close-up view of the ratio using the expression with the 
first three terms in (39). This clearly demonstrates the remarkable precision of the resurgence relation (39): 
the agreement is better than the one percent level already by fc « 10. 


we deduce the numerical value of the Stokes constant Si: 

Si Ri —0.040883 f , for AA = 4 parameters (40) 

Note that Si is pure imaginary. 

^ The factor 27ri is the denominator is conventional, being convenient for the associated Borel analysis. 
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The form of the large order growth is intimately connected with the Borel transform of the 
asymptotic series The Borel transform is defined as 


k=0 


AO) 

Jk+1 k 

k\ ^ 


(41) 


For factorially divergent series such as this new function has a finite radius of 

convergence. The Borel summation of {w) is then formally defined as 

(w) = + / e-"'* (s) ds . (42) 

Jo 

The angle 9 determines the contour of integration in the complex- s plane, usually called the “Borel 
plane”, and is correlated with the phase of the original expansion parameter w. For our original 
problem we are interested in 0 = 0, since w is a real parameter. However has singularities 

along the real axis, where 9 = 0. A quick way to see this is to use the large order expression (39) 
in (41): 


/^°^(s)~^r(/3 + i)(5-s)-^-i 


f(i) I HI) I (^-«)^ HI) I 

JO Q Jo o! o T\Jo 


+ 


(43) 


/? ' /3(/3-l)' 

It is clear that there is a branch point at s = S'. In order to avoid the branch cut we deform the 
contour infinitesimally as 0 = 0^, but depending on the direction, one encounters a discontinuity. 
Let us consider 9 = 0^. The discontinuity is pure imaginary, and is of the order Im[f^^\w)] oc 
Si/q^^ Requiring that the trans-series is real, this ambiguous contribution to that 

arises from the singularity in the Borel transform must be canceled by This cancellation 

fixes the imaginary part of a in the trans-series (30) to be [12, 13] 


Im[a] 


-^Xm[Si] 


(44) 


With J\f = 4 parameters, Im{Si) ~ —0.04. For a detailed analysis of the reality of trans-series, 
and the associated cancellations of imaginary ambiguities we refer the reader to [13]. 

The branch cut singularity of the Borel transform f^^\s) can also be seen in another way. By 
generating many coefficients fj^\ using the recursion relations (18), we automatically generate 
many terms in the Taylor series for from (41). The location of the nearest singularity 

may be deduced by a root test or ratio test, but to see the branch cut it is better to use a Bade 
approximant [33]. The Bade approximant to f^^\s) approximates it as a ratio of two polynomials: 


/<“>(») 


p(0)(5) 

g(‘^)(s) 


(45) 
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The branch cuts of the actual Borel transform function manifest themselves as an accumulation of 
poles in the Fade approximant. Of course the further away a branch point is from the origin, the 
more precision and number of terms it would require to reproduce the associated branch cut. We 
computed the symmetric Fade approximant of order 300; namely (s) and (s) are taken to 
be polynomials of order 150. We used precision of 800 significant figures. The result is shown in 
Fig. 2. It is clear that there is a branch cut that starts at s = 5* = along the positive real 
axis, consistent with the resurgent behavior in (43) and the large-order behavior in (39). The fact 
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FIG. 2: The poles in the Borel complex s-plane of the Fade approximant to the Borel transform of the 
hydrodynamic expansion, The poles accumulate into a branch cut that starts at S' = 5 ^) which is 

the analog of the instanton action. 


that the location of the singularity in the Borel plane coincides with the “instanton action” (31) 
in the trans series is a generic feature of resurgence. 

Due to the nonlinear nature of the original differential equation (15), there are in fact infinitely 
many exponentially suppressed terms in the trans-series, each of which is associated with an action 
nS. This fact translates into the existence of branch points located at s = nS for all n > 1. 
The ambiguities that arise from each of these branch points are cured by the existence of the 
non-hydrodynamic series. In terms of the large order growth this means that the late terms of the 
hydrodynamic series actually contains information about all the non-hydrodynamic series. More 
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precisely, the expression (39) can be further refined as [13] 


AO) 

Jk 




Si 

+s 

+s 


nk + f3) / ( 1 ) 

27ri5fc+/3 yo 

2 r(fc + 2/3) / 

^ 2TTi (25)^+/^ \ 

3 r(fc + 3/3) / 
^ 27ri (35)^+^ \ 


fc + Zl-l-'i (fc + /3-l)(fc + /3-2) ^2 

(2) 2S (2) _(25f_ (2) 

A: + 2/3-l-'i (/s + 2/3-l)(A; + 2/3-2) ■'2 

(3) 35 (3) _(35f_ (3) 

A: + 3/3-l-'i (/s + 3/3-l)(A; + 3/3-2) ■'2 



+ ... . 


( 46 ) 


Note that the expression in each line involves low order coefficients of non-hydrodynamic series 
of different order, and that there is only one constant, the Stokes constant Si, that needs to be 
determined numerically. 


B. First (leading) non-hydrodynamic expansion: 


We can repeat the analysis of the previous section for the fluctuations {w) about the first 
non-hydrodynamic term in the trans-series. Once again, general arguments for a real trans-series 
predict the large order growth [13]: 



~ 25i 


r(fc + /i) 

2Tri 


( fO) _|_ ^ _ fO) _ ^ _ d 2 ) 

V° k + /3-l-'^ {k + /3-l)ik + l3-2)-'^ 


(47) 


( 2 ) 

The coefficients ^ on the right-hand side are the low order terms of the second non-hydrodynamic 
series f^^\w) in (36). Note that in (47) all constants on the right-hand-side are known: the overall 
normalization constant is hxed to be twice the very same Stokes constant Si found in the large- 
order behavior (39) of the hydrodynamic series. In this sense (47) is an even stronger prediction 
than (39). In Figure 3 we plot the ratio of the large order expression (47) to the exact coefficients 
generated from the recursion relations (33). Again, the agreement is excellent. 

The Borel plane structure can be deduced by the same Borel-Pade analysis, and as shown in 
Figure 4 we see once again a cut along the positive real axis, starting at the branch point S. The 
Borel plane structure is almost identical to that of the hydrodynamic series; namely there are 
infinitely many branch cuts located at s = nS and accordingly, the large order growth (47) can be 
refined similar to (46). In other words, late terms of contain the low order terms of all the 

higher non-hydrodynamic series /^”^(rc) with re > 2. 
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FIG. 3: The first plot shows plot the ratio of the large order expression (47) to the exact coefficients 
as a function of k. The blue, green and red points refer to the inclusion of the zeroth, first, and second 
subleading terms in (47). The second plot shows a close-up view of the ratio using the expression with the 
first three terms in (47). This clearly demonstrates the remarkable precision of the resurgence relation (47): 
the agreement is at the one percent level already by fc « 10. 
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FIG. 4: The poles in the Borel complex s-plane of the Fade approximant to the Borel transform of the first 
non- hydrodynamic expansion, /0)(s). The poles accumulate into a branch cut that starts at S' = 2 ^- 

C. Second (next-to-leading) non-hydrodynamic expansion: 

At the next order, studying the large order growth of the coefficients of the second non- 
hydrodynamic series we observe a new phenomenon, the possibility of which was pointed 

out by Aniceto and Schiappa in their exhaustive analysis of the resurgent structure of real trans- 
series [13]. The large order growth of the ' coefficients is determined not only by the the low 
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order coefficents of but also by the the low order coefhcents of 

i-Sf 


f(2) 

Jk 


Tik-/3) / (1) -5 (1) 

27ri (-5)^-/^ ^ ^ 


{k-/3-l){k-l3-2) 




+ 


3Si I /f > + 


27ri 


5 


A; + /3- 1 


/f ^ + 


^2 


(k + /3-l)(k + /3-2) 


/f+ ... 


+ . 


( 48 ) 


Here is a new Stokes constant whose value we determine numerically to be S'-! ~ —57.922 + 
115.651z. Note that in the first line the phase of S-i is cancelled by the factor i(—l)~^ in the 
denominator, as the coefficients are real. 

In Fig 5 we plot the ratio of the expression (48) to the exact coefficients generated from 
the recursion relations in (35). As in the previous cases, the agreement is excellent. Note that 
if we do not include the effect of the coefficients, the agreement is terrible. In fact, for the 
particular choice of AA = 4 parameters that we study, because /? < 0, the contribution of the first 
non-hydrodynamic sector (i.e. the first line in (48)) is more dominant compared to the third non- 
hydrodynamical sector (the second line in (48)). Related with this fact, the coefficients of the series 
are actually sign-alternating, as opposed to the hydrodynamic series and the first 

non-hydrodynamic series for which the expansion coefficients are non-alternating. But this 
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FIG. 5: The first plot shows plot the ratio of the large order expression (48) to the exact coefficients , 
as a function of k. The blue, green and red points refer to the inclusion of the zeroth, first, and second 
subleading terms in (48). The second plot shows a close-up view of the ratio using the expression with the 
first three terms in (48). This clearly demonstrates the remarkable precision of the resurgence relation (48): 
the agreement is at the one percent level already by fc « 10. 

sign-alternating behavior does not necessarily mean that is Borel summable. To explore 

this, we use the Borel-Pade method to study the singularity structure of the Borel transform (s). 
In the complex-s Borel plane (see Fig. 6 ), in addition to the infinitely many branch points along the 
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positive real axis s = nS, associated with the higher (3'’'^, 4*^, 5*^, etc.) non-hydrodynamic series, 
there is an additional branch cut along the negative real axis that starts at the branch point s = —S. 
This branch cut is associated with the first non-hydrodynamic series, and the corresponding Stokes 
constant is related to the Stokes discontinuity at 0 = vr. It is clear from the existence of the 



FIG. 6: The poles in the Borel complex s-plane of the Fade approximant to the Borel transform of the 
second non-hydrodynamic expansion, /(^^(s). The poles accumulate into two branch cuts, which start at 
S = 2 ^, and at S' = 

branch cuts on the positive real axis, is not Borel summable along the positive real axis, 

even though it is an alternating series. A better way to conceptualize this is to view 
as being composed of different sub-series, each of which is related to different non-hydrodynamic 
sectors (each line in (48)). The sub-series related to the first non-hydrodynamic series is alternating, 
Borel summable along the positive real axis, and is the dominant contribution to the late term 
coefficients Therefore the whole series is alternating. However, all the sub-series 

associated with non-hydrodynamic series with n > 3 lead to non-alternating sub-series which are 
not Borel summable along the positive real axis. This phenomenon is reminiscent of the quantum 
mechanical example discussed in [34], where there are complex instantons with negative actions, 
dubbed “ghost instantons”, which contribute to the large-order behavior of perturbation theory; 
the structure is very similar to the way in which the first non-hydrodynamic series contributes to 
the large order growth of the ' expansion coefficients. 
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D. Third and higher non-hydrodynamic expansions: with n > 3 


The large-order behavior found in the second non-hydrodynamic sector in the previous section 
persists at third order and beyond. For example, at third order we find the large-order growht of 
the expansion coefficients: 


f(3) 

Jk 


2S-1 


45i 


r(fc-/?) 

27ri 

nk±R( M) 


/f+ 






{k-^-l){k-l3-2) 


_ rw j_ 

2mS^+Py^ k + P-l 


S 


fy + 


52 


{k + /3-l)ik + /3-2) 


# + •■• ) + 




(49) 


with the same S-i parameter. The agreement of (49) with the exact coefficients generated from 
the recursion relations (35) is again excellent, as shown in Figure 7. We stress that there are no free 
parameters in this comparison. All constants on the right-hand-side of (49) are determined. The 
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FIG. 7: The first plot shows plot the ratio of the large order expression (49) to the exact coefficients , 
as a function of k. The blue, green and red points refer to the inclusion of the zeroth, first, and second 
subleading terms in (49). The second plot shows a close-up view of the ratio using the expression with the 
first three terms in (49). This clearly demonstrates the remarkable precision of the resurgence relation (49): 
the agreement is at the few-percent level already by /c « 10. 


Borel plane structure of the Borel transform is shown in Figure 8, showing the same two 

branches, starting at s = ± 2 ^' expect that this Borel plane structure persists to all orders. 

For completeness, we note that there are further even smaller corrections to (49), involving 
and also with n > 5. The part introduces a new Stokes constant S -2 which can be 
determined numerically. The general pattern for the large order growth of the hydrodynamical 
series is as follows. There is a set of alternating sub-series related to the sectors n = 

1,2,..., N — 1, each of which has an associated Stokes constant 5_(jv-n)' They generate branch 
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FIG. 8: The poles in the Borel complex-s plane of the Fade approximant to the Borel transform of the third 
non-hydrodynamic expansion, The poles accumulate into a two branch cuts that start at s = 

and s = — 2 ^- The two poles that do not lie on the real axis are artifacts of the finite order (150) of our 
Fade approximation. 

points in the Borel plane along the negative real axis: s = —S, —25,..., —{N — 1)5. As mentioned 
before, there is another set of sub-series related to the sectors n = iV -|- 1, Ai -|- 2,... all of which 
are non-alternating. They generate branch cuts at s = 5, 25,.... In other words the location of 
the branch cut in the Borel plane of (s) due to the non-hydrodynamical sector n ^ N is given 
by the “relative action” Sbranch = (n — N)S. This is an illustration of the general behavior of real 
trans-series and has also been observed in resurgent analysis of topological string theory [7, 13]. 


VI. TRANS-ASYMPTOTIC REARRANGEMENT AND INITIAL CONDITIONS 


At earlier proper-time (i.e. smaller w), the terms in the trans-series expansion (29) become 
disordered, as exponentials compete in size with inverse powers of w. In such a situation, there 
is a systematic way to rearrange the terms of the trans-series as one approaches the edge of its 
domain of numerical usefulness [1, 35]. This is the first step in “trans-asymptotic matching” [35]. 
This procedure effectively resums all the exponentially small non-hydrodynamic modes for a given 
power of 1/w, and rearranges the trans-series (30) into the form 


f{w) ~ 

k=0 


wk 


(50) 
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This means that we have identified Fk{cr Q with the formal expansion 


Fk{(TO = ^ 


( 51 ) 


71=0 


but we will see that we can in fact obtain closed-form expressions for the Fk{<jC,). 

In (50), C = C,{w) is the same “instanton fugacity”, C('f^) = introduced in (31). 

Inserting this reorganized ansatz for f{w) into the original differential equation (15), leads to 
a sequence of ordinary differential equations, one for each Fk{C), with (( now being regarded as 
the independent variable. Note that the argument of Fk involves also an arbitrary numerical 
factor multiplying (^. As in the expansion (30), there is just one such undetermined constant, and 
by comparison we identify it with the trans-series parameter a in (30). The tower of differential 
equations for Fk{C) can be solved recursively: the first equation of the tower involves only To(C), the 
second only involves To(C) and Fi{Q, and so on. [For notational convenience we scale a = 1 here, 
and reintroduce it again later]. Moreover, while the equation for Fq{Q is nonlinear, all subsequent 
equations for Fk{C) are linear. Furthermore, when expanded as series in (], each function Fk{C) is 
convergent, even though the full trans-series expression (50) is of course still divergent [35]. 

It is straightforward to show that the ansatz (50) leads to the following equations for Fk{C)- 

Cr E (sc^ + - ik' - 1)T,7_i(C) ) Fk-k'iO + E Fk'iOFk-k'iO 

k'=0 \ i.,_n 




k'=0 

k-l 


l)Fk-i{C)) + ^CrFk-iiC) + 4a E Fk-k'-iiOFk'iC) 

^ k'=0 

+ (l + FkiO + g(a - C,)4-1,0 + 3 (l + <^^>0 = 0 (52) 

We focus on the class of parameters with Cx = 2C^, which includes the particular case of AA = 4 
parameters that we have been using for numerical purposes. Then the first equation 


- ^Ci"o(C)^^ - IcFoiC) + HFoiOf + 3Fo(C) + ^ = 0 


(53) 


determines Fo(C) to be 

Fo(C) = ^(-2 + cV±cCV4+^) 




(54) 


where c is an integration constant. Matching with the trans-series expansion (29) fixes c = 3 and 
selects the upper sign: 

3 


1 


i"o(C) = -x + xC + CA I + tC 


9 


(55) 
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From the small C expansion of -Fo(C) can read off the leading (i.e. w^) coefficient of the 
non-hydrodynamic series. To see this explicitly, compare 


= /(S°^ + C + ^ C' + + (56) 

with the leading coefficients given in Eqs. (19), (34), (36) and (37), with Cx = 2(7,,. 

The next equation in the tower (52) determines Ti(C) 

~ {P '0 + 1) + 3 (2Fo -|- 1) Ti — 3Crf C (To -|-1) -|- ACtFq + -CtFq + g (C't ~ C*)?) = 0 

(57) 


Note that this equation is linear in Fi{(^). After choosing the integration constant to match the 
1/w term of the first non-hydrodynamical expansion (34), we obtain 

( 3C + V4 + 9C2 : ' 

T’i(C) = 


(3C+y/4T^ 


54(7„(a - C’„)C + Cr (9((7r - Cr^)Q^ + 2C„) (-3C + 

i(X ^ 


+QCT{2Cr — 3Cr]) ( sinh 


(58) 


Similarly to (56), the small C expansion of Fi{(^) generates all the coefficients in the trans-series: 


i'.s; 9^3 


10 - 


9C, 

a 


C + 13(7,, — 


18(72 


a 


+ 4a r 




= + /f ^ c + /f ^ c' + /f ^ c" + 


(59) 


which can be checked from Eqs. (19), (34), (36) and (37), with Cx = 2(7,,. This procedure can be 
continued to determine all the Fk{C), as all the equations are linear for A; > 1. 

The advantage of this procedure can be seen when we reintroduce the trans-series parameter 
a factor in the argument of F^, writing again Fk{crQ. The condition of reality of f{w) fixes the 
imaginary part of a in terms of the (numerically determined) Stokes parameter Si in equation (40, 
44), but the real part of a is arbitrary [13, 35]. This remaining real constant parametrizes the entire 
family of trans-series solutions to the differential equation, each member of which connects to the 
hydrodynamical expansion at large w, but which have different behavior once w becomes small 
enough that the exponential factors become comparable with powers of w. Thus, this real part of 
the trans-series parameter characterizes the various different possible “initial conditions”, at early 
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f(w) 



FIG. 9: The black (dashed) line is the numerical integration of the nonlinear differential equation (15), 
starting at very early time w = 10“®. The red (solid) lines correspond to the trans-asymptotic expansion 
(50), using the expressions for Fq and Fi in (55, 58), the imaginary part of the trans-series parameter a 
fixed by (40, 44), and different values of the real part of a (from bottom to top: —0.2, —0.12,0.06,0.12, 0.2). 
This clearly illustrates that the physical meaning of the real part of the trans-series parameter a is that it 
corresponds to matching a particular member of the family of trans-series to a particular solution with a 
prescribed initial conditon. A detailed example of such a matching procedure for the Painleve I equation 
can be found in [36]. 

times (i.e. small w), each of which eventually tends to the hydrodynamical expansion as tc —)• oo. 
In [15] this phenomenon was referred to as an “attractor”. This is in fact a generic feature of the 
large class of nonlinear differential equations whose solutions can be expressed in trans-series form 
[1, 35]. This behavior is illustrated in Figure 9, in which we vary the real part of cr, and see that 
this corresponds to different possible choices of initial condition. This simple observation resolves 
the apparent inconsistency between the fact that the formal hydrodynamic series expansion (19) 
has no undetermined parameter, while the exact (numerical) solution to the first-order nonlinear 
equation (15) clearly depends on one initial condition: this is why there is a family of trans-series 
solutions, parametrized by the trans-series parameter a, more precisely by its real part if f{w) is 
real. 


25 






VII. CONCLUSIONS 


We have shown that the nonlinear ordinary differential equation describing boost-invariant 
Bjorken flow in conformal hydrodynamics exhibits characteristic signs of resurgence. The formal 
late-time Taylor expansion of the hydrodynamic derivative expansion is asymptotic, and can be sys¬ 
tematically extended to a trans-series expansion, for which we have characterized the Borel plane 
structure. This analysis reveals precise relations between the fluctuations about different “non- 
perturbative sectors”, associated with modes beyond hydrodynamics. These modes are exponen¬ 
tially damped and the damping rate, or the inverse relaxation time, plays the role of an instanton 
action in the language of semi-classical physics. The factorial large order growth of the fluctuation 
coefficients in one sector is precisely determined by low orders of expansions about neighboring 
sectors. This means, in particular, that the asymptotic hydrodynamic expansion encodes physical 
information about non-hydrodynamic modes. For example, the trans-series expansion of the en¬ 
ergy density £{t) at late proper-time includes also exponentially small terms ~ r“ exp [—cr^/^], 
and these non-hydrodynamic modes have an interpretation in terms of quasi-normal modes in the 
gravitational picture. Resurgent asymptotic analysis suggests that these general features should 
extend to the partial differential equations in more general hydrodynamical problems with gravity 
duals, in which the metric coefficient functions are expanded in terms of both proper time and 
distance from the horizon. 
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